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ANALYTICAL APPROXIMATION FOR STEADY SHIP WAVES 
AT LOW FROUDE NUMBERS 


Francis Noblesse 
David Taylor Naval Ship R&D Center 
Bethesda, MD 20084 


ABSTRACT 

A simple analytical relationship between a ship- 
hull form and its steady far-field Kelvin wake is obtained 
by considering the low-Froude-number limit of the 
Neumann-Kelvin theory. In particular, this relationship 
predicts the occurrence of a sharp peak in the amplitude 
of the waves in the far-field Kelvin wake at an angle, a, 
from the ship track that is smaller than the Kelvin-cusp 
angle of 19°1/2 for a hull form which has a small region 
of flare and is wall sided elsewhere, if the Froude 
number is sufficiently small. An explicit relationship 
between the angle, @, between the ship track and the 
tangent to the ship mean waterline in the region of flare 
and the corresponding ‘‘wave-peak’’ angle a in the 
Kelvin wake is obtained. For instance, this relationship 
predicts the occurrence of a sharp peak in wave 
amplitude at an angle a in the Kelvin wake equal to 14° 
for a hull having a small region of flare within which the 
waterline-tangent angle @ is approximately equal to either 
30° or 74°. This theoretical result may explain the bright 
returns that have sometimes been observed in SAR 
images of ship wakes at angles smaller than the Kelvin- 
cusp angle. The low-Froude-number asymptotic analysis 
of the Neumann-Kelvin theory presented in this study 
also predicts that the wave-resistance coefficient is O(F2), 
where F is the Froude number, for a ship form with a 
region of flare, O(F*) for a ship form that is wall sided 
everywhere and has either a bow or a stern (or both) that 
is neither cusped nor round, and O(F®) for a wall-sided 
ship form with both bow and stern that are either cusped 


or round. 


1. INTRODUCTION 

It has been observed, see for instance Fu and Holt 
(1982) and McDonough et. al. (1985), that SAR 
(Synthetic Aperture Radar) images of ship wakes 
sometimes reveal bright returns along rays at angles from 
the track of the ship smaller than the Kelvin-cusp angle 
of 19°1/2. A plausible explanation for these surprising 
observations was proposed by Scragg (1983) who 
considered a simple ship bow form with large flare for 
which he found that the zeroth-order slender-ship 
approximation to the far-field wave-amplitude function 
given in Noblesse (1983) predicted a sharp peak in the 
value of the amplitude of the divergent waves at an angle 
from the track of the ship equal to approximately half 
the bow entrance angle. This numerical result of Scragg 
was confirmed by Barnell and Noblesse (1986) who also 
found that the peak in the amplitude of the divergent 
waves becomes sharper as the value of the Froude 
number decreases, and thus suggested that the occurrence 
of a sharp peak in the amplitude of the far-field Kelvin 
waves was a large-flare low-Froude-number feature. 

The numerical studies of Scragg and of Barnell 
and Noblesse are based on two simple approximations to 
the far-field wave-amplitude function, namely the Michell 
thin-ship approximation for which no peak was found 
and the zeroth-order slender-ship approximation which 
exhibited a peak as was already noted, so that it is not 
clear from these studies whether a more realistic 
mathematical model for the far-field wave-amplitude 
function, such as that provided by the Neumann-Kelvin 


theory, would also predict the occurrence of peaks in the 


amplitude of the far-field Kelvin waves. Furthermore, the 
numerical results obtained by Scragg and by Barnell and 
Noblesse correspond to a particular ship form and thus 
provide little physical insight into the origin of the 
predicted peak in the amplitude of the far-field Kelvin 
waves, specifically the manner in which such a peak is 
related to the shape of the ship hull. 

A complementary analytical study of the low- 
Froude-number limit of the Neumann-Kelvin theory for 
an arbitrary ship form is thus presented here. This 
asymptotic analysis of the Neumann-Kelvin theory 
provides a simple analytical relationship between a ship- 
hull form and its steady far-field Kelvin wake. In 
particular, this-relationship predicts the occurrence of a 
sharp peak in the amplitude of the waves in the far-field 
Kelvin wake at an angle, a, from the ship track that is 
smaller than the Kelvin-cusp angle of 19°1/2 for a hull 
form which has a small region of flare and is wall sided 
elsewhere, if the value of the Froude number is 
sufficiently small. A simple explicit relationship between 
the angle, w, between the ship track and the tangent to 
the ship mean waterline in the region of flare and the 
corresponding wave-peak angle a in the Kelvin wake is 
given and depicted in Figure 3b. For instance, this figure 
predicts the occurrence of a sharp peak in the amplitude 
of the divergent or transverse waves at an angle a in the 
Kelvin wake equal to 14° for a hull having a small region 
of flare within which the waterline-tangent angle @ is 
approximately equal to 30° or 74°, respectively. This 
analytical result may explain the bright returns that have 
sometimes been observed in SAR images of ship wakes 
at angles smaller than the Kelvin-cusp angle. 

The low-Froude-number asymptotic analysis of the 
Neumann-Kelvin theory presented in this study also 
predicts that the wave-resistance coefficient is O(F?), 
where F is the Froude number, for a ship form with a 
region of flare, O(F4) for a ship form that is wall sided 
everywhere and has either a bow or a stern (or both) that 
are neither cusped nor round, and O(F) for a wall-sided 
ship form with both bow and stern that are either cusped 
or round. 

More precisely, the low-Froude-number asymptotic 
approximation (50) to the far-field wave-amplitude 
function K(t) shows that the main contributions to the 


function K(t) stem from several particular points on the 


mean waterline. These are the bow and the stern, on one 
hand, and (usually but not always) one (or several) 
point(s) of stationary phase. Indeed, the number of these 
points of stationary phase, and their position on the 
waterline, depend on the value of t and on the shape of 
the waterline. The first two terms in the low-Froude- 
number asymptotic expansions for the contributions Kgs 
of the bow and stern and the contributions K + of the 
points of stationary phase in equation (50) are given by 
equations (51)-(55) and (59)-(62), respectively. The 
second-order terms in these asymptotic expansions are 
defined by complex expressions. However, the first-order 
terms provide simple approximations defined explicitly in 
terms of the geometrical characteristics of the hull and 
the velocity components in the tangential directions t and 
NX t to the hull (see Figure 2). In particular, these low- 
Froude-number asymptotic expansions show that the 
contributions K, and Kg of the bow and stern are O(1) 
except if the bow or stern is cusped or round, in which 
case we have Kgs = O(F?). The contribution of a given 
point of stationary phase is O(1/F), and thus is 
dominant, if the hull has flare at this point; otherwise, 
that is if the hull is wall sided at the point of stationary 
phase, its contribution is O(F). The low-Froude-number 
approximation (50) also shows that we have K(t) = 

O(1/t3) as t > ©. In fact, this result is valige? any 


value of the Froude number. 


The latter result implies that the lines along which 
the steepness of the short divergent waves in the far-field 
Kelvin wake takes given large values, say 1/7 and 1/15, 
are parallel to the ship track, as was found in Figure 21 
of Barnell and Noblesse (1986) by using the Michell thin- 
ship approximation for a simple ship form. The 
Neumann-Kelvin theory therefore predicts that the far- 
field Kelvin wake contains three distinct regions: (i) a 
narrow constant-width inner region bordering the track 
of the ship where no divergent gravity waves can exist, 
(ii) an outer region where the usual transverse and 
divergent waves are present, and (iii) an intermediate 
region at the boundary between the inner and outer 
regions where short steep divergent waves can be found. 
In reality, surface tension and possibly also viscosity 
must evidently be taken into account in the vicinity of 
the track of the ship. 


2. STATIONARY-PHASE APPROXIMATION TO 
THE FAR-FIELD KELVIN WAKE 


The far-field Kelvin wake may be conveniently 
analyzed in terms of the nondimensional far-field 
coordinates (7,6) = (X,Y,Z) g/U2, where g is the 
gravitational acceleration, U is the speed of advance of 
the ship, and (X,Y,Z) are dimensional coordinates. The 
mean free surface is taken as the plane € = 0, with the 
¢ axis pointing upwards, and the & axis is chosen along 
the track of the ship, that is in the ship centerplane, and 
pointing towards the bow. The origin of the system of 
coordinates is placed within the ship. 

Equations (4), (3a), (7) and (8) in Barnell and 
Noblesse (1986) yield the following expression for the 
elevation ¢ of the free surface at a sufficiently large 
distance behind the ship, such that nonlinearities may be 


neglected: 


n€(En) = Re di 


where K(t) is the far-field wave-amplitude function and 


% (E, +E_) K(t) (1+t?)!/2 dt, (1) 


E,, is the exponential function defined as 


E,(t;£,0) = exp [i£6,(t;0)], (2) 
with the phase-function 6, and the parameter o defined 
as 

6,(t;0) = (1Fot)(1 + t?)!/2, (3) 
o = n/(-&). (4) 


For a ship with port- and starboard-symmetry, as is 
considered here, the Kelvin wake is symmetric about the 
ship track n = 0. We may then restrict the analysis of 
the Kelvin wake to the quadrant n 2 0, — < 0 and 
assume o 2 O. Let a be the angle between the track of 
the ship and the line joining the origin of the system of 
coordinates to the observation point (€,n). We thus have 


a = tan~'o. (5) 


A far-field asymptotic approximation, valid as — > 
—o, to the wave integral (1) may be obtained by 
applying the method of stationary phase, as is well 
known. The result of this classical asymptotic analysis 
may be found in Barnell and Noblesse (1986), for 
instance. Specifically, equations (28), (24a), (20b), (25a-d) 
and (26a-d) in this reference yield 


(—&)!/7t(E.a) ~ Re 
{A_K(t_) exp [i(E@_ —1/4)] 


+ A,K(t,) exp [i(EO, +1/4)]} as + —™, (6) 


with 0 << < tan710—3/2) m 19°28", 


where t,(a) are the values of t for which the phase- 
function 8, (t;a) is stationary, ©,(a) = 6, (t,:) 
represents the corresponding values of the phase-function 
8, (t;a), and A +(a) is the function defined as A (a) = 
(2/m)'/2(1 + t,,2)!/2/[£0'/(t,:0)]'/2. The functions t,(a), 
©,(a) and A (a) are given by 


t,(a) = [1+(1—80%)!/2}/4o, (7) 
@,(a) = [3(1 —80%)!/2} 

[1 + 40+(1 — 807)!/2)!/2/21/2g5, (8) 
A,(a) = [1 +407+(1 — 802)!/293/421/47 


4n!/293/2(1 ~~ 802)!/4, (9) 
where we have 
o = tana (10) 


as is given by equation (5). 

The far-field asymptotic approximation (6) shows 
that the wave pattern at any point (,a), with << —1 
and 0 < a < 19°28’, consists in two elementary plane 
progressive waves, so-called transverse and divergent 
waves, as is most well known. The wavelengths A, and 
the directions of propagation B,, measured from the 
track of the ship, of the transverse (A_,f_) and 
divergent (A, ,8,) waves at an angle a from the track of 
the ship are given by 


A(a) = 2!/216n07/[3=(1 — 80%)!/7] 


[1 — 407+(1 — 807)!/2)!/2, (11) 
B.(a) = cos !{2!/220/[1 + 407+(1 — 807)!/2]!/2}. (12) 
We have 
1 >A_/2n > 2/3 >A, /2n > 0 and (13) 
0<p_ <sin~'(1/3!/?) ~ 35°16’ < B, < 90°. (14) 


Equation (6) shows that the amplitudes of the 
transverse and divergent waves in the Kelvin wake are 
asymptotically equal to A, |K(t,)|/(—£)!/* as + —&. 
The steepnesses, say s,, of these waves then are given by 
s, = A,|K(t,)|/A,(—£)'/2. We then have 
(-8)!s (Ea) ~ S,|K(t,)| as § > —=, (15) 
where the functions S,(a) are defined as S,(a) = 
A,(a)/A, (a). Equations (9) and (11) then yield 
S,(a) = [3(1 — 807)!/4] [1 — 407+(1 — 807)!/2)!/2 


[1 + 402+(1 e 802)!/2]3/4 721/46 4n3/207/2(1 by 8021/4, (16) 


Equations (7), (8), (9), (10) and (16) show that we 


have 

t_= 0,0. = 1, A_= (2/m)!”2 (17a,b,c) 
and S_ = 1/2!/2n3/2 for a = 0, (17d) 
t, © 1/2a, 0, ~ 1/4a, A, 1/2n!/2q3/2 (18a,b,c) 
and S,~ 1/16n/2a7/? as a > 0, (18d) 
t, = 1/2! and ©, = 3(3/2)!/274 (19a,b) 
for a = tan~ '(1/23/2) ~ 19°28’, 

A, © 3/n!/261/4(1 — 802)!/4 and (19c) 
S, ~~ 9/413/261/4(1 — 802)!/4 (19d) 


as a > tan !(1/23/2). 


The stationary-phase values t,(a), the phase-functions 
8_,(a), the amplitude-functions A +(a) and the steepness- 
functions S,(a) are depicted in Figures 1a and 1b, which 
correspond to the transverse and divergent waves, 
respectively. 

It may be seen from these figures and from 
equations (19c) and (19d) that the amplitude-functions 
A,,(a) and the steepness-functions S + (a) are singular at 
the Kelvin cusp line o* = 1/8, @ = tan '(1/23/2) 
19°28’, in accordance with the well-known fact that the 
asymptotic approximation (6) is not uniformly valid at 
the boundary of the Kelvin wake. A complementary 
asymptotic approximation, expressed in terms of Airy 
functions, valid at and near the Kelvin cusp line is given 
in Ursell (1960) for the particular case of the Kelvin wave 
pattern due to a concentrated pressure point at the free 
surface. However, Ursell’s more complex asymptotic 
approximation will not be considered here because the 
simple asymptotic approximation (6) is little affected by 
the weak singularity (1 —802)~ !/4 for points (E,a) inside 
the Kelvin wake and not too near the cusp line a 
19°28’, in which we are mostly interested in this study. 

Equation (11) yields 


A, © Bna? as a > 0. (20) 


This approximation and the approximation (18b) show 
that there are an infinite number of divergent waves with 
indefinitely shorter wavelength in the vicinity of the track 
of the ship, as is well known. It may be seen from 
Figure 1b and equations (18a,c,d) that the stationary- 
phase value t , (a), the amplitude-function A , (a) and the 
steepness-function S , (a) are unbounded in the limit a > 
0. Equations (18a,c,d) yield 


Oe 3s? 6 9 12° 15°" 182 


Qa —— > 


Fig. la - The Stationary-Phase Value t_(a), the Phase- 
Function @_(a), the Amplitude-Function A_(a) and 
the Steepness-Function S_ (a) Corresponding to the 

Transverse Waves in the Kelvin Wake 
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Fig. 1b - The Stationary-Phase Value t , (a), the Phase- 
Function © , (a), the Amplitude- -Function A , (a) and 
the Steepness-Function S,(@) Corresponding to the 

Divergent Waves in the Kelvin Wake 


A, © (2/n)!/t , 3/2 and (21a) 
S, vt, 7/2/2!/2n/2 ast, = ©. (21b) 


Equations (21a) and (6) show that the amplitude of the 
divergent waves, given by A , |K(t , )|/(—£)!/”, vanishes 
at the track of the ship if we have 


t3/2K(t) > Oast > ~. (22) 


Furthermore, it is shown in Barnell and Noblesse (1986) 
that the asymptotic approximation (6) is uniformly valid 
at the track of the ship if condition (22) is verified. It 
may be seen from equations (21b) and (15) that the 
steepness of the divergent waves, given by 

S, |K(t,)|(—2)!”, is unbounded at the track of the 
ship if we have 


t’/2|K(t)| > © ast > ~. (23) 


Both conditions (22) and (23) can be satisfied 


simultaneously if we have 
|K(t)| ~ 1/t# as t > © with 3/2 <p < 7/2. (24) 


In summary, the asymptotic approximation (6) 
expresses the far-field wave pattern of a ship at a point 
(E.a), with £ < —1 and 0 < a < 19°28’, as the sum of a 
transverse wave and a divergent wave. The phase 
£0 _(a)+n/4 of these two waves are defined explicitly in 
terms of — and a; and their amplitudes 
A,(a)K(t,)/(-8)!/? are given by the product of the 
functions A,(a)/(—£)!/*, which are also defined 
explicitly in terms of é and a, and the far-field wave- 
amplitude function K[t,(a)], which depends on the speed 
(Froude number) and the shape of the ship. The far-field 


wave-amplitude function is now considered. 


3. BASIC EXPRESSIONS FOR THE FAR-FIELD 
WAVE-AMPLITUDE FUNCTION 


The far-field wave-amplitude function K(t) may be 
conveniently defined in terms of the nondimensional 
near-field coordinates (x,y,z) = (X,Y,Z)/L, where L is 
the length of the ship. In the Neumann-Kelvin theory, 
the function K(t) is given by the sum of an integral 
around the mean waterline of the ship and an integral 
over the mean wetted-hull surface. Specifically, for a 
ship with port and starboard symmetry, as is considered 


here, the function K(t) may be expressed in the form 

K(t) = K,(t) + K_(), (25) 
where the functions K, are given by 

F°K,(t) = le E,(n,7+1,4,—n,t,¢4 + ivpé)t, dl 


+ fi, exp(v?p*z)E,(n, + vp*$n,) da, (26) 


as is shown in Noblesse (1983). In this equation, F is the 
Froude number defined as 


F = U/L)”, (27) 
v is its inverse, that is 

m = IVA, (28) 
p is defined as 

p = (1+t?)!”?, (29) 
and E, represents the exponential function 

E, = exp [—iv*p(x+ty). (30) 


Furthermore, c and h represent the positive halves of the 
mean waterline and of the mean wetted-hull surface, 
respectively. The unit vector tangent to c and pointing 
towards the bow is denoted by t(t,,t,,0), and Tin, n,n) 
is the unit vector normal to h and pointing into the 


xrlys 
water, as is indicated in Figure 2. The term n, is defined 
as 


n, = —n,+i(n,+tn,)/p. (31) 


Figure 2 - Definition Sketch for a Single-Hull Ship with 
Port and Starboard Symmetry 


Also, dl and da represent the differential elements of arc 
length of c and of area of h, respectively. Finally, 
>= (x) represents the nondimensional disturbance 
potential ¢ = ®/UL at the integration point x on c or h, 
$, represents the derivative of $ in the direction of the 
tangent vector t to c, and $4 is the derivative of > in the 
direction of the vector nxt, which is tangent to h and 
pointing downwards as is shown in Figure 2. ° 
Equations (25) and (26) express the far-field wave- 
amplitude function K(t) in terms of the value of the 
Froude number and the form of the mean wetted-hull 
surface, as was noted previously. More precisely, 
equation (26) expresses the function K_(t) as the sum of 
a line integral around the mean waterline and a surface 
integral over the mean wetted-hull surface of the ship. 
Furthermore, these integrals involve the disturbance 
velocity potential > in their integrands. The relationship 


between the far-field wave pattern of a ship and its form 


and speed thus is a complex one, which offers little 
physical insight. However, analytical approximations for 
the waterline and the hull integrals in equation (26) can 
be obtained in the limiting case when the value of the 
Froude number is sufficiently small. This low-Froude- 
number asymptotic approximation for the function K(t) 


is now obtained. 


4. LOW-FROUDE-NUMBER LAPLACE 
APPROXIMATION TO THE FAR-FIELD WAVE- 
AMPLITUDE FUNCTION IN TERMS OF A - 
WATERLINE INTEGRAL 


For large values of v = 1/F, or more generally of 
vp, the exponential function exp(v2p2z) in the hull- 
surface integral in equation (26) vanishes rapidly for 
negative values of z. Therefore, only the upper part of 
the mean-wetted hull surface h yields a significant 
contribution in the low-Froude-number limit. More 
precisely, the hull-surface integral can be approximated 
by a line integral around the mean waterline c. 
Furthermore, this integral can be combined with the 
waterline integral in equation (26). The analysis is 
presented in detail in Noblesse (1986a) and is briefly 
summarized below. 

The mean waterline is represented by the 


parametric equations 
X = Xq(A) and y = y((A), (32a,b) 


where the parameter A varies between its bow and stern 


values, that is 
Ap SA < Asc. (32c) 


In the vicinity of the mean free surface, the hull surface 


72) 


represented by the parametric equations 


X = X9(A)+z x,(A)+27x,(A)+..., (33a) 
y = yolA)+z y,(A)+z7y,A)+..., (33b) 
where Ap <A < Ag and z <0. (33c) 


The velocity potential $(A,z) on the hull surface in the 
vicinity of the plane z = 0 likewise is expressed in the 


form 
$ = $o(A) +z $0) +2790) +... (34) 


Differentiation of the functions x,(A), y,(A), ¢,(A) with 
respect to the parameter A is denoted by the 
superscript '; thus, we have xj = dxg(A)/da. 

By using the foregoing parametric representations 


for the hull surface and the velocity potential, applying 


the Laplace method for approximating the hull-surface 
integral in equation (26), and combining the resulting 
waterline integral with that already present in equation 
(26), we may obtain ____ after lengthy algebraic 
transformations ____ the following low-Froude-number 
approximation for the function K(t): 


K(t) ef SE *a, +E) a_)dAasF-+0, (35) 
where we ee 
q = I/p = 1/(1 + t?)!72, (36) 
E,~ is the exponential function 
E,* = exp[—ivp(xptty,)], (37) 
and a, is the amplitude function 
a, = u,a,*+Fq*(u,)a,* + O(F4), (38) 
where u, is defined as 
u, = I/[I —iq(x,+ty,)], (39) 
and the functions a,* and a,~ are now defined. 
The first-order amplitude function a,* is given by 

ay* = yo’ A,/(1 +62) + 2q(X9'+tyo’)(u,)B, $y 

+ C$’ + uD,$,/(1 +e?) + ip(y,4o’ —Yo'$,), (40) 
where we have 
€ = (Yo X,—Xq y,)/u, (41) 
with u defined as 
UW = [%/ + Oo) alls (42) 


and the coefficients A,, B,, C, and D, are defined as 


A, = [(1+pyp'/u)(1 —pyp'/u) +e?) 
+ i(pyg’ /u)[y, (xp +tyo')/u te}, (43a) 
B, = Q(yyFtx,) + i(yyX_—X,Y>), (43b) 
Cc, = [Gt- P2X9'Yo' /u?) 
+ (Pyp'/u)€(x, Xp" +yY1¥q')/u(1 + &’)] 
+ ily,(X9'+tyg’)/u+e] 
[PXq'/u— (Pq! /U)e(X,Xq' +¥1¥q')/u( +e)],  (43c) 
D, = [(%p/#tyg')/ull(1 + €?)(y, +iat)u, + i(pyo’ /wey,] 
— (PYg'/u)e(Pyg’ /u — ie). (43d) 
The second-order amplitude function a,* is given 
by 
ay = —y,' + 2u,(dom,* +4 )m,* + $2Mo* —i¥q'¥2”) 
+ 6i(u,)7[¢9(m,*y.* + My~¥37) ats $)M9¥271 


12(u,)*4omo*(y27), (44) 


where we have 


Yn = Q(x,+ty,), (45) 
m,* = u, + iq(y,'Ftx,'), (46) 
with 

Bopecls (47a) 
Hy = XyYy' YX" + 2EQY9' —Y2%X'), (47b) 


Ha = Xo! —YyXq" + 2%)" —yoxXy') 


+ 3(X,Y9' —¥3Xp')- (47c) 


In the particular case when the phase of the 
exponential function La is stationary, that is at a point 
(X9,¥o,0) where we have Xo'ttyo’ = 0, the first-order 
amplitude function a,~ takes the form 


ay = + eliyg’ + (py, *it)$o' Fud,]/(1+ie) 
+ IP(¥j#' —¥o'4)) if Xp'+tyy’ = 0. (48) 


Furthermore, if the hull surface intersects the plane z=0 


orthogonally at a point of stationary phase we have 
ay = i F’pyy’ 874,8x? 
if Xo +tyo’ = 0 and i, = 0, (49) 


and the amplitude function a, then is or order F2. 

The low-Froude-number asymptotic approximation 
to the far-field wave-amplitude function K(t) given by the 
waterline integral (35) is considerably simpler than the 
exact expression (26), which involves both a waterline 
integral and a hull-surface integral; and the approximate 
expression (35) is well suited for efficient numerical 
evaluation. However, expression (35) can be simplified by 
applying the method of stationary phase, which takes 
advantage of the rapid oscillations of the exponential 
function E,~ given by equation (37) in the low-Froude- 
number limit v > ©, or more generally in the limit 
vp > ©. This stationary-phase approximation is now 


obtained. 


5. LOW-FROUDE-NUMBER STATIONARY-PHASE 
APPROXIMATION TO THE FAR-FIELD WAVE- 
AMPLITUDE FUNCTION 


The method of stationary phase indicates that the 
major contributions to the integral (35) in the limit when 
the exponential functions E,~ are rapidly oscillating, that 
is if v > © or more generally vp ~ ©, stem from points 
where the phases of these exponential functions are 
stationary, that is from points where we have X,'+tyo’ = 


0, and from the end points Ap and Ax of the integration 


range, that is from the bow and the stern. This 
stationary-phase analysis is presented in detail in 
Noblesse (1986b) and only its results are given here. We 
have 


K() ~ iq(K,-K, += K,) as F — 0, (50) 


where K, and K, correspond to the contributions of the 
bow and the stern, respectively, and K, corresponds to 
the contribution of a point where the phase of the 
exponential function En is stationary, that is where we 
have yo’ /(—X9') = +1/t or dyp/dXxp = ¥1/t; the 
summation in the asymptotic approximation (50) is 
extended over all such points of stationary phase on the 
mean waterline c. The expressions for the contributions 
of the bow and stern Kgs and of the points of 
stationary phase K_, are given below. 

The contribution Kgs of the bow or stern may be 


expressed in the form 
Kgs = Ag sexp(— iv*pxp 5), (51) 


where Xp ¢ is the abscissa of the bow or stern and the 


amplitude-function Ags Is given by 
Ags = 2n,A, +iF’q(A,* +A,~)+O(F4); (52) 


in this expression, the functions A, and A, are defined 


by the equations 


(n, + iqn,)A, =n,n, [1 —ipn,(n, + iqn,)/(1 — rite — p*n,”)] 


+ te - ipnytyeq(1 - n,? +ipn,n,)/(1 — n, = p*n,’) (53) 
and 
@,.'A,* = (u,a,*/0,')’ +igtu,)2a,*, (54) 


where the superscript ‘ denotes differentiation with 
respect to the parameter A, the functions u,, ay and a,* 
are defined by equations (39), (40) and (44), respectively, 


and the function 6, is defined as 


6, = Xp+tyg. (55) 


+ 


Expression (54) for the second-order amplitude- 
function A, is a complex one. However, equation (53) 
defines the first-order amplitude function A, explicitly in 
terms of the value of t, p = (1 +t2)!/2 and q = 1/p, the 
geometrical characteristics of the hull at the bow or 
stern, namely the unit vector t(t,,t,,0) tangent to the 
mean waterline and the unit vector Tn (n,,ny,n,) normal 
to the hull, and the components +, and $4 of the velocity 
vector in the directions of the unit vectors t and nxt 
tangent to the hull. Equation (52) shows that the first- 


order approximation to the amplitude-function Ags 


given by 2n,A,, vanishes if n,= 0, that is if the 
waterline has a cusp at the bow or stern. Equation (53) 
shows that the first-order approximation to the function 
Ags also vanishes if the bow or stern is round, since we 
then have Wt = O and ¢, = =A. = 0 by symmetry. It 
may thus be seen that the contribution of the bow or 
stern to the far-field wave-amplitude function K(t) is of 


order F? if n, = Oorn, = 0 at the bow or stern, 


y 
respectively; that is, we have 


Kgs = O(F’) ifn, = 0 or ny = 0 (56) 


at the bow or stern. In the particular case when the hull 
surface is vertical at the bow or stern we have n, = 0 
and equation (53) becomes 


A, = t,—$,—ipt,¢,/(I—pt,’) ifn, = 9, (57) 
which yields 
A, = t,-$,+O(F*) ifn, = 0 (58) 
as is indicated by the free-surface boundary condition , 
= -F4,, = OF). 

The contribution K., of a point of the mean 


waterline c where the phase of the exponential function 
Eo* is stationary, that is where we have 


dyp/dxy = ¥1/t, (59) 
may be expressed in the form 
K, = +v(2nr)!/2A, 

exp [Fiv7(yot, — Xpty)/ty + ien/4], (60) 
where r is the radius of curvature of c at the point of 
stationary phase (Xp,yo), € is equal to +1 or —1 if the 
center of curvature of c at the point (Xx,Yo) is upstream 


or downstream from (Xp,Y), respectively, and the 


amplitude-function A, is given by 
A, = (1—n,2)!/2(n,t, + $4) — F2aA,* + O(F4); (61) 
in this expression, the second-order amplitude-function 
A,~ is defined by the equation 
+2[(x9')? + (¥9')7]!/27A5* =a [aya y/Oxcali 

+u,a,*[5(0, /"")?/30,"" — 6,4)/4(,"")? 

+ 2iq(u,)?a,*, (62) 
where the superscript ‘ denotes differentiation with 
respect to the parameter A and the functions u,, a,*, a,* 


and 6, are given by equations (39), (40), (44) and (55), 
respectively. 


The expression for the second-order amplitude- 
function A,* is a complex one. However, equation (60) 
and the first-order approximation to the amplitude 
function, namely 


A, = (1=n,2)!2(n,1, + $4) + OF), (63) 


provide a simple explicit expression for the stationary- 
phase contribution K, in terms of the geometric 
characteristics of the hull and the downward tangential 
derivative ¢, of the potential at the point of stationary 
phase. In the particular case when the hull surface is 
vertical at the point of stationary phase we have n,=0 
and equation (63) becomes 


A, = —$,+O(F’) ifn, = 0, (64) 
which yields 
A, = O(F*) ifn, = 0. (65) 


Equation (60) then shows that the stationary-phase 
contribution at a point (X,Y) Of c where the hull is 


vertical is of order F,. that is we have 
K, = O(F) ifn, = 0. (66) 


On the other hand, equations (60) and (63) show that we 
have 


K, = O(1/F) ifn, #0. (67) 


The stationary-phase contribution at a point where the 
hull has flare thus is dominant in the zero-Froude- 
number limit. 

The summation in equation (50) is extended to all 
the points of the mean waterline c where the phase of 
the exponential function E)* or the function Ep ~ is 
stationary, that is the points where the slope dyo/dXxp of 
c is equal to —1/t or +1/t, respectively. The number of 
stationary points, and their position along the waterline, 
depend on the value of t and on the shape of c. For 
instance, for the simple case of a hull with waterline 
consisting of a sharp-ended parabolic bow region 1/4 < 
x < 1/2 defined by the equation y = 4bx(1 —2x), where 
b denotes the ship’s beam/length ratio, a straight parallel 
midbody region — 1/4 < x < 1/4, and a round-ended 
elliptic stern region — 1/2 < x < —1/4 defined by the 
equation y = b[—2x(1 +2x)]!/2, there is one point of 
stationary phase in the stern region given by x = 
—[1+1/(1 + 4b*t2)!/2]/4, so that we have —1/2 <x < 
—1/4 for 0 <t < © with x > —1/2 ast > 0 andx > 


—1/4 as t ~ ©, and one point of stationary phase in the 


bow region given by x = (1+1/4bt)/4 for 1/4b < t < 
oo, so that we have 1/2 > x 3 1/4 for 1/4b € t € &© 
with x > 1/2 ast > 1/4b and x > 1/4 ast > ©. We 
thus have one point of stationary phase in the stern 
region for 0 < t < 1/4b and two points of stationary 
phase, one in the stern region and one in the bow region, 
for 1/4b < t < ~. The two points of stationary phase 
approach the shoulders x = +1/4, where dy/dx = 0, as 
li == Coe 

The asymptotic approximation (50) and equations 
(51)-(58) and (59)-(67) defining the contributions of the 
bow and stern and of the stationary-phase point(s) on 
the waterline, respectively, show that the low-Froude- 
number behavior of the far-field wave-amplitude 
function is strongly influenced by the shape of the hull in 
the vicinity of the waterline. More precisely, for a value 
of t for which there is one (or more) point of stationary 
phase on the mean waterline where the hull has flare, the 
contribution of this stationary-phase point dominates the 
contribution of the bow and stern:and is of order 1/F, 
that is we have K(t) = O(1/F). On the other hand, for a 
value of t for which either there corresponds no point of 
stationary phase or the hull has no flare at the point(s) 
of stationary phase, the dominant contribution stems 
from the bow and stern, and it is of order 1, that is we 
then have K(t) = O(1). However, if n, = 0 or ny = 0 
at both the bow and the stern, that is if the bow and the 
stern are either cusped or round, their contribution is 
O(F2) and the contribution of the stationary-phase 
point(s), which is O(F) if there is no flare (as is assumed 
here), is dominant; so that we then have K(t) = O(F). 

For a ship form that is everywhere wall sided, the 
contribution of the bow and stern is dominant for all 
values of t, and we have K(t) = O(1) forO0 <t < ~. On 
the other hand, for a hull form that has flare over a 
portion of the waterline and is wall sided elsewhere, the 
contribution of the bow and stern is dominant, and O(1), 
only for those values of t for which the corresponding 
points of stationary phase fall outside the range of flare; 
for the range of values of t for which the corresponding 
points of stationary phase are within the range of flare, 
the contribution of these stationary-phase points is 
dominant, of order 1/F. In this instance, the function 
K(t) is O(1/F) for a range of values of t (corresponding 
to the region of flare) and O(1) for other values to t. If 
the region of flare is of small extent, and the slope 
dy/dx of the waterline does not vary widely within that 


region, the range of values of t for which K(t) = O(1/F) 
is also small. The far-field wave-amplitude function K(t) 
can then exhibit a sharp peak for some value of t in the 
low-Froude-number limit. In fact, several isolated peaks 
of the function K(t) can exist if the hull form has several 
distinct regions of flare within which the slope of the 
waterline varies gradually. 

It should be noted that the result K(t) = O(1/F) 
for values of t for which the hull has flare at the 
corresponding points of stationary phase does not imply 
that the corresponding free-surface elevation becomes 
unbounded in the zero-Froude-number limit F = 0. 
Indeed, the asymptotic approximation (6), where we have 
(E,0)=(X,Z)g/U%, then yields Zg/U2 = O[1/F(—£)!/2] = 
O[1/(— X/L)!/2] as F > 0 and —X/L — ©. The free- 
surface elevation Z thus is of order (U2/g)/(— X/L)!”2 as 
F > 0 and —X/L > ©, and equation (15) shows that 
the corresponding wave steepness s is O[1/(— X/L)!/2}. 
For values of t for which K(t) = O(1) it is seen that 
Zg/U? and s are O[F/(— X/L)!/2] as F > 0 and — X/L 
een OOs 

It should also be noted that the asymptotic 
approximation (51)-(55) for the contribution of the bow 
and stern is not uniformly valid for the values of t for 
which the bow or the stern is a point of stationary 
phase, that is for which the waterline slope dyp/dxp at 
the bow or the stern is equal to —1/t or 1/t, 
respectively. Indeed, it may be shown that we have 
1—n,?—p’n,? = 0 at a point of stationary phase, so 
that the first-order approximation to the amplitude 
functions Ags given by equation (53) becomes 
unbounded. The asymptotic approximation (60)-(62) for 
the contribution of a point of stationary phase likewise is 
not uniformly valid at a stationary-phase point where th 
waterline has an inflexion point. Indeed, the radius of 
curvature r at such an inflexion point is infinite and 
equation (60) yields an unbounded contribution Kor 
Asymptotic approximations valid for these special cases 
may easily be obtained and are given in Noblesse 
(1986b). It will only be noted here that the far-field 
wave-amplitude function K(t) for a hull form having 
flare at a point where the waterline has an inflexion may 
be expected to exhibit a particularly pronounced peak at 
the value of t corresponding to the inflexion point since 
we have K(t) = O(1/F4/3) as F > 0 for this particular 


value of t. 


6. CONCLUSION: HULL FORM AND KELVIN- 
WAKE FEATURES 


The classical far-field asymptotic approximation to 
the Kelvin wake, obtained in section 2 by applying the 
method of stationary phase, and the low-Froude-number 
asymptotic approximation to the far-field wave-amplitude 
function, obtained in sections 4 and 5 by successively 
using the Laplace method and the method of stationary 
phase, provide a simple analytical relationship between 
the hull shape, on one hand, and the waves it generates, 
on the other hand. This explicit relationship between the 
wavemaker and its waves is summarized below. 

The far-field asymptotic approximation to the 
Kelvin wake (6) shows that at any point (£,a), with & = 
Xg/U2 « —1 and 0 < a < tan~!(2~ 3/2) w 19°28’, the 
wave field consists in two plane progressive waves ___ a 
transverse wave and a divergent wave _____ with 
wavelengths A_ and A, propagating at angles #_ and 
6B, from the track of the ship, respectively. The 
wavelengths A, and the propagation angles f., depend on 
the angle from the ship track a alone, that is A, and p, 
are independent of the hull shape and size, as is well 
known. Specifically, the functions A(a) and fi,(a) are 
defined by equations (10), (11) and (12). At a given 
downstream distance £, the amplitudes of these waves, 
on the other hand, are given by the product of the 
functions A,(a), defined by equation (9), and the far- 
field wave-amplitude function K(t) evaluated at the 
stationary values t(a) given by equation (7). The 
function K(t) depends on the hull shape and the Froude 
snumber in a fairly complicated manner via an integral 
over the mean wetted-hull surface and an integral around 
the mean waterline, as is indicated by equation (26). A 
low-Froude-number asymptotic approximation to these 
integrals is obtained in sections 4 and 5. 

The analytical approximation (50) shows that for a 
given value of t corresponding to a given value of a, as 
is specified by equation (7), the main contributions to the 
function K(t) stem from several particular points on the 
mean waterline. These are the bow and the stern, on one 
hand, and (usually but not always) one (or several) 
point(s) of stationary phase. Indeed, the number of these 
points of stationary phase and their position on the 
waterline, defined by the condition 


|dyy/dxy| = 1/t, (68) 


depend on the shape of the waterline and the value of t. 
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Let @ denote the angle between the tangent to the 
mean waterline and the track of the ship, that is we have 


tan p = |dyg/dxp| and 0 < @ < n/2. (69) 


Equations (68), (69), (7) and (10) then show that the 
stationary values t , (2) and t_(qa) associated with a given 
value of a are defined by the equation 


@,(a) = tan~!{40/[1+(1 — 80?)!/2}}, (70) 


where o = tana. For a given waterline shape, equation 
(70) thus defines the number of stationary points and 
their position on the waterline corresponding to any 
given angle a inside the Kelvin wake. In particular, 


equation (70) yields 
O<, <tan !(2!/2) ~ 54°44’ <p < 90° (71) 


and m, © 2a, p_ ~ n/2-aasa— 0. (72a,b) 


Points of the waterline with slope between 0 and 54°44’ 
thus contribute mostly to the system of divergent waves 
while waterline slopes between 54°44’ and 90° mostly 
contribute to the transverse waves. 

Equations (10), (11), (12) and (70) define the wave- 
lengths ,, the wave-propagation angles /3, and the 
waterline-tangent angles m, corresponding to a given 
angle from the ship track a. The functions A,(a), B,(a) 
and @,(a) are depicted in Figure 3a, where the subscripts 
T and D are used, instead of — and +, to refer to the 
transverse and divergent waves, respectively. The 
foregoing relationships between a and A, f, @ may be 
used for determining the angle from the ship track a, the 
wavelength A and the wave-propagation angle # 
corresponding to a given waterline-tangent angle q. 
Specifically, we may obtain the remarkably simple 


relations 
a = tan! [tanp/(2 +tan2)], (73) 


A/2n = sin’ and B = n/2-@. (74a,b) 


The functions a(q), A(p) and B(@) are depicted in Figure 
3b. Alternatively, the foregoing relationships among 4, a, 
i and f# can be represented in the form of Figures 3c and 
3d, which depict the functions (jf), a(f3), A(B) and (A), 
a(A), B(A), respectively. These equivalent graphical 


representations show that we have 


0 < @p < tan” (21/2) w 54°44" < , < 90°, (75a) 
0 < Ap/2n < 2/3 < Ay/2n < 1, (75b) 
90° > Bp > tan !(2~ 1/7) w 35°16' > fy > 0, (75c) 
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Fig. 3a - The Waterline-Tangent Angle @, the Wave- 
Propagation Angle B and the Wavelength A as Functions 
of the Angle From the Ship Track a; the Subscripts T 
and D Refer to the Transverse and Divergent Waves, 
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Fig. 3c - The Waterline-Tangent Angle », the Angle 
From the Ship Track a and the Wavelength A as 
Functions of the Wave-Propagation Angle p. 


where the subscripts D and T refer to the divergent and 
transverse waves, respectively. 

The first two terms in the low-Froude-number 
asymptotic expansions for the contributions Kgs of the 
bow and stern and the contributions K,, of the points of 
stationary phase in equation (50) are given by equations 
(51)-(55) and (59)-(62), respectively. The second-order 
terms in these asymptotic expansions are defined by 
complex expressions. However, the first-order terms 
provide simple approximations defined explicitly in terms 
of the geometrical characteristics of the hull and the 
velocity components in the tangential directions t and 
nxt to the hull. In particular, the low-Froude-number 
asymptotic expansions given in section 5 show that the 
contributions Ky and K, of the bow and stern are O(1) 


except if the bow or stern is cusped or round, in which 
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Fig. 3b = The Angle From the Ship Track a, the Wave- 
Propagation Angle B and the Wavelength A as Functions 
of the Waterline-Tangent Angle q. 
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Fig. 3d - The Waterline-Tangent Angle @, the Angle 
From the Ship Track a and the Wave-Propagation Angle 
B as functions of the Wavelength A. 


case we have Kgs = O(F2). The contribution of a given 
point of stationary phase is O(1/F), and thus is 
dominant, if the hull has flare at this point; otherwise, 
that is if the hull is wall sided at the point of stationary 
phase, its contribution is O(F). 

Thus, for a ship form that is everywhere wall 
sided, the contributions of the bow and stern are 
dominant (assuming that they are not both either round 
or cusped) for all values of the angle a from the ship 
track, that is everywhere in the far-field Kelvin wake. On 
the other hand, for a hull that has flare over a portion 
of the waterline and is wall sided elsewhere, the 
contributions of the bow and stern are dominant, and 
O(1), only for those angles a in the Kelvin wake for 
which the corresponding points of stationary phase on 


the waterline fall outside the region of flare; for the 


range of values of a for which the corresponding points 
of stationary phase are within the range of flare, the 
contribution of these stationary-phase points is 
dominant, of order 1/F. In this instance, the amplitude 
of the waves in the far-field Kelvin wake is of order 
(U2/g)/(— X/L)!/2, as F > 0 and X/L > —©, for the 
range of values of a corresponding to the region of flare 
and of order F(U2/g)/(— X/L)!/2 for values of a outside 
this range. If the region of flare is of small extent, and 
the waterline-tangent angle @ does not vary widely within 
that region, the corresponding range of values of the 
angle a where the wave amplitude is an order of 
magnitude larger than elsewhere is also small, and thus 
appears as a peak for sufficiently small values of the 
Froude number. This peak is particularly pronounced for 
a hull with a small region of flare in the vicinity of an 
inflexion point of the waterline. 

The low-Froude-number asymptotic analysis of the 
Neumann-Kelvin theory presented in this study thus 
shows that the characteristics of the far-field Kelvin wake 
strongly depend on the shape of the ship hull, notably 
the presence of flare and the shape of the waterline at 
the bow and stern. This analysis also predicts that the 
nondimensional wave-resistance coefficient R/eU2L?, 
where U and L are the speed and the length of the ship 
and @ is the density of water, which is given by the 


Havelock integral 


foe} 
mR/oU2L2 = F4 ifs \K(t)|2(1 +t?) !/2 dt, (76) 


is O(F2) for a ship form with a region of flare, O(F*) for 
a ship form that is wall sided everywhere and has either 
a bow or a stern (or both) that is neither cusped nor 
round, and O(F®) for a wall-sided ship form with both 
bow and stern that are either cusped or round. 

The low-Froude-number asymptotic approximation 


(50) also shows that we have 


K(t) = O(1/t3) as t > ©. (77) 


This result is actually valid for any value of the Froude 
number; indeed, the asymptotic approximation (50) is 
valid not only in the low-Froude-number limit F > 0 but 
more generally in the limit vp = (1 +12)!/2/F — oo, that 
is as F ~ 0 or/and as t > ©, as may be seen from tne 
exponential functions exp (v2p2z) and E, = 
exp [— ivp(x+ty)] in equation (26). 

Equations (15), (18a), (21b) and (77) then yield 


(SYA S Ga) Yh NEO WOES => 0 (78) 
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for the steepness of the short divergent waves in the 
vicinity of the track of the ship. By using equations (4) 
and (5), which yield a ~ n/(—&) as a > 0, in equation 
(78) we may obtain 


s, (En) © 1/n!? as n > 0. (79) 


Equation (79) thus shows that the lines along which the 
steepness of the short divergent waves in the far-fieid 
Kelvin wake takes given large values, say s, = 1/7 and 
1/15, are parallel to the ship track, as was found in 
Figure 21 of Barnell and Noblesse (1986) by using the 
Michell thin-ship approximation for a simple ship form. 
The Neumann-Kelvin theory therefore predicts that the 
far-field Kelvin wake contains three distinct regions: (i) a 
narrow constant-width inner region bordering the track 
of the ship where no divergent gravity waves can exist, 
(ii) an Outer region where the usual transverse and 
divergent waves are present, and (ili) an intermediate 
region at the boundary between the inner and outer 
regions where short steep divergent waves can be found. 
In reality, surface tension must evidently be taken into 


account in the vicinity of the track of the ship. 
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